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Abstract 

We show that the /o(980) and ao(980) states appearing in scattering and production processes 
can be described quite well by the two-channel Oller-Oset-Pelaez (OOP) version of unitarized chiral 
theories. It is impossible, however, to deduce a decisive conclusion on the nature of them from the 
fitting to experimental data alone. Using explicit N c dependence of parameters in OOP amplitudes 
we demonstrate that light scalar resonances are not of qq states but dynamical rescattering effects 
generated under chiral symmetry, unitarity and channel couplings. 

@ @ @ @ @ @@ @ @ 

1 Introduction 

Probably we have not yet reached a consensus on the nature of light scalar mesons below 1 GcV[l, 2]. 
Main issues around them are whether they are real propagating particles or are dynamical rescattering 
effects, whether they are of qq states, four quark states or KK molecule states, and whether they form a 
nonet within themselves or only one or two are needed to form a different nonet with higher mass scalar 
states. 

In order to study these issues many authors focus on the analyses of the /o(980) and ao(980) states 
appearing in decay processes such as <f> — > 7/0 (980) and 4> — > 700 (980). In most theoretical approaches the 
radiative <p decays take place through a radiative charged kaon loop diagram so as to satisfy the Okubo- 
Zweig-Iizuka (OZI) rule. The difference in various approaches comes from the way how to construct 
amplitudes K + K~ — > rri\mi including the /o(980) or ao(980) states, where m\ and m-i are final two 
pseudoscalar mesons. Among the ways we quote models using / (980) and ao(980) propagators [3, 4, 5], 
unitarized chiral approaches[6, 7, 8, 9, 10], a model using the linear sigma modelfll, 12], a compound 
model of chiral 0(p 2 ) terms and propagators[13], and a model with channel couplings between hadronic 
and quark channels [14]. There are approaches without using the explicit kaon-loop mechanism; a model in 
which production amplitudes are assumed to be proportional to scattering amplitudes with real coefficient 
functions in order to satisfy unitarity [15], and a quark model with a radiative quark loop by assuming 
that all these mesons are to be dominantly qq states [16]. Similarly, 77 — > f and ao processes are 
expected to reveal the quark structure of them [17] . How unambiguously can we determine both theoretical 
and experimental branching fraction B(<f> — > 7/0) ? Calculating a pole and its residues of the /o(980) 
state, Boglione and Pennington [15] have given B(<j) — > 7/0) much smaller than the usually quoted value 
(3.4 ± 0.4) x 10" 4 [18]. This indicates that we should be cautious of extracting the correct signal of the 
/o(980) state from experimental mass distributions. If this is the case, it becomes obscure whether we 
can reach a decisive conclusion on the nature of the /o(980) and ao(980) states through these analyses. 

In this paper we attempt to describe the /o(980) and ao(980) states appearing in scattering and 
production processes within the two-channel Oller-Oset-Pelaez version[19], that is an approximate version 
of Inverse Amplitude Method (IAM)[20, 21, 22, 23, 24] to unitarize amplitudes up to 0(p 4 ) of Chiral 
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Perturbation Theory (ChPT)[25, 26]. The basic ingredients of SU(3)xSU(3) ChPT are the octet Nambu- 
Goldstone particles ir, K, rj s , where we treat r]$ as the physical rj in the two-channel OOP version, and any 
scalar meson fields are not introduced as independent degrees of freedom in advance. It offers, therefore, 
a good theoretical framework to study the issues. Of course, this does not imply that the existence 
or non-existence of a resonance can be predicted by IAM, since ChPT amplitudes of 0(p 4 ) include a 
set of phenomenological parameters, low energy constants (LECs) , which are to be determined so as to 
reproduce experimental data. After having fixed the LECs, the multi-channel amplitudes are determined 
without any additional free parameters. Though our set of the LECs are not the best solution, the phase 
shifts are reproduced quite well. Since two-meson mass spectra of the radiative </> meson decays depend 
on the off-diagonal K + K~ — > irir and K + K~ — > ttt] scattering amplitudes below the KK threshold, the 
resulting mass spectra offer a proving ground of the OOP amplitudes. It turns out that the mass spectra 
are reproduced fairly well by the OOP amplitudes. In 77 — ► tttt processes, however, the / (980) signal 
is masked by complicated interferences, since various production mechanism contribute coherently. This 
process is a proving ground of production mechanisms, if we have more accurate data near the /o state, 
therefore. 

Even if IAM or its OOP version can give simultaneous fits to scattering and production data, it 
is impossible to deduce a decisive conclusion on the quark structure of the /o(980) and do (980) states 
because of the reason stated above on the LECs. We should note, however, that the pion decay constant 
/,r and the LECs in ChPT have a definite N c dependence, where N c is the number of colors. Due to 
this nature of the theory we can study how scalar amplitudes behave when N c increases from 3. It has 
been shown recently that the scalar states do not survive as N c is larger than 3, and then the scalar 
states are probably not qq states but dynamical ones, while vector mesons survive as narrow resonances 
as typical qq states[27, 28]. In this sense IAM and its approximate OOP version can give not only a 
unified description of scattering and production processes but also an insight into the quark structure of 
the light scalar mesons. 

We discuss the OOP description of scattering processes in the next section, including a discussion on 
LECs used in this paper. Scalar meson productions in radiative 4> meson decays and 77 collisions are 
discussed in Sec. 3, and the large N c behavior of scalar states are discussed in Sec. 4, and concluding 
remarks are given in the last section. 

2 The OOP description of scattering processes 
2.1 The OOP amplitude 

The ingredients of IAM consist of chiral order 0(p 2 ) and 0(p 4 ) amplitudes of ChPT. An 0(p 2 ) amplitude, 
denoted by T^ 2 ^(s, t, u), has a form of a linear function of s, t, u and meson mass squared divided by f 2 , 
where s, t and u are the usual Mandelstam variables. An 0(p 4 ) amplitude is composed of polynomial 
terms and one-loop terms: A polynomial term, denoted by T^ ly (s,t,u), is given as a sum of quadratic 
functions of s, t, u and meson mass squared with the LECs, denoted follows: 

T$ y (s,t,u)= J2 jiL n P n (s,t,u), (2.1) 

n=l,8 

where P n are the quadratic functions. Loop terms are given by x with divergent loop integrals, 
which are regularized in the MS — 1 scheme at the renormalization scale /J.[25, 26]. An s-channel loop 
term is given by t^- 2 \s)J(s)t^ 2 \s), where is a partial wave amplitude derived from T^(s,t,u) and 
J(s) is the one-loop function which will be given soon. There are similar t- and w-channel loop terms and 
tadpole terms, but they are discarded in the OOP version. The ingredients of amplitudes in the OOP 
version consist of the 0(p 2 ) terms, the polynomial and the s-channel loop terms of 0(p 4 ) amplitudes. 
We use the 0(p A ) amplitudes given in Ref. [23], where the kaon and r\ decay constants are set equal to 
the pion decay constant /„• in order to guarantee exact perturbative unitarity. 

The two-channel OOP amplitudes are written in the symmetric 2x2 matrix form as 

T(s) = t (2 \s)[t^( S ) - tW(8)]-HW(s), (2.2) 
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where t^ 2 \t^) is the 2x2 partial wave amplitude with chiral order 0{p 2 ) (0(p 4 )) amplitude and s is 
the total center of mass (CM) energy squared. Our T-matrix is normalized as 



S^s) = 5 ii -2ipt'\s)T ij {s)p]'\s), (2.3) 

1 9h 

"« = ^Ts d(S - Si) i2A) 

where pi is the phase space factor with ki (sj) being the CM momentum (threshold energy squared) of 
the z-th channel. The partial wave (s) is given as 

t^( S )=C y ( S )+t^( S )-J( S )-t^( S ), (2.5) 

where (s) is the partial waves of (s, t, u), and J(s) is the diagonal 2x2 loop integral; Ji(s) of 
the i-th channel is written 



1 f / mim 2 \ A 2 i / m 2 

-— j < -1 + log — + log 

{Any [ v m / s V m i 



Vl-(miim 2 ) 2 /s, (2.7) 
Ai(«) = o-i+(*)ai_(a) (2.8) 

for the channel with unequal masses mi and m 2 in the i-th channel, and A 2 i = to 2 ~ m i- Because the 
imaginary part of Ji(s) is given as 

ImJi(s) = -pi(s), (2.9) 
T(s) satisfies the exact s-channel unitarity relation 

ImT^s) - -2£(«)p fc (s)71y(s). (2.10) 

The LECs, L\ to is, are determined so as to reproduce S'-wave tttt phase shifts up to about 1.2 GeV. 
It should be noted that our fitting covers a wide energy range from the tttt threshold to 1 GeV or more 
including resonances, that the LECs appear in T non-linearly and that the OOP version discards the t- 
and u-channel loop terms and tadpole terms of the full 0(p 4 ) amplitudes. So that our set of the LECs and 
p could be different from those of ChPT determined through low energy data below the lowest resonance 
region. The scale p has a meaning corresponding to a cut-off parameter as discussed in Ref. [19], and we 
set p = 0.87 GeV, that is larger than usually adopted value p = m p by about 100 MeV. Our set of the 
LECs are given as 



L 1 = 0.70 x 10~ 3 , L 2 = 1.30 x 10~ 3 , L 3 = -3.20 x 10" 3 , 
L 5 = 1.50 x 10~ 3 , L 7 = -0.25 x 10~ 3 , L H = 0.71 x 10~ 3 ' ; ; 



and L4 and L e are fixed equal to zero. 



2.2 The (J, J) = (0, 0) channel: mr x KK 

The characteristic behavior of this channel is that the 7T7t phase shift <5j° rises from the tttt threshold, 
forms a plateau of 50° - 80° from 500 to 800 _MeV, crosses 90° in the region 800 - 900 MeV, and 
suddenly increases over 200° just below the KK threshold. This feature is reproduced quite well as 
shown in Fig. 1 (a). The phase shift calculated in the single tttt channel formalism cannot exceed 90° 
as shown by the dotted line in Fig. 1 (a). On the other hand a bound state pole appears in the KK 
channel, when the channel coupling to the tttt channel is switched off. A weak coupling to the tttt channel 
generates a narrow resonance behavior in the tttt channel. The steep rising of the tttt phase shift near the 
KK threshold is due to the interference between the 7T7T background and the narrow resonance born as 
the KK bound state. In this sense the /o(980) state is called a bound state resonance[29, 30, 31, 32, 33]. 
Also shown are the production rate defined by (1 — ?7q )/4 with 7700 being the inelasticity and the the 
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phase, 5®2 — ^11 + ^22 1 above the K K threshold. These give the size and the phase of the off-diagonal 
scattering amplitude T12 above the KK threshold. 

In our set of the LECs and \i the phase shift Sf® crosses 90° at 875 MeV. The scattering lengths are 
given as 

a 00 = 0.217 m- 1 a 20 = -0.041 m" 1 . (2.12) 
According to Ref. [42] ChPT up to two loop accuracy gives 

a oa = 0.220 ±0.005 a 20 = -0.0444 ± 0.0010. (2.13) 




Figure 1: (a) The irir phase shift S^. The solid (dotted) line represents the result of the two- (single 
7T7r-)channel calculation. The vertical (horizontal) line represents the KK threshold (90°). Experimental 
data are full circles [34], blank squares [35], bars [36, 37], and triangles [38]. (b) irir — > KK phase, S^. 
Experimental data are taken from [39]. (c) (1 — r/Q )/4. The data are full circles[39] and bars[40]. (d) <5 20 
with the data from [41]. The abscissas are the invariant tttt mass in units of GeV. 



The phase difference c>( 00 ) — 5( 20 ) at the kaon mass is 49.9°, which is the phase of two-pion decay amplitudes 
of a kaon originating from the final state interactions. Recently the KLOE collaboration gives (48±3)° [43] , 
and ChPT does (47.7 ± 1.5)° [42]. 




Figure 2: Energy dependence of —Tij. (a) —Tu. (b) — T12. (c) — T22. The solid lines are the absolute 
values, the dotted lines the real parts, and the dashed lines the imaginary parts. 



We also show all of the amplitudes in Fig. 2. The /o structure appears as a sharp dip in Tu 
owing to the interference with the large background, which is often called the cr(600) state, but T 12 and 
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T 22 show the clear peak corresponding to the /o state just below the KK threshold. The magnitude and 
shape of Ti 2 below the KK threshold is crucial to the <j) 1™ decays, and it turns out that our OOP 
amplitude T i2 reproduces the nir mass spectra fairly well. 

2.3 The (I, J) = (1, 0) channel: nr) x KK 

Since there are almost no experimental data on extrapolated tttj elastic scattering, the ttt] mass distri- 
butions from K~p — > 7r _ ^S+(1385) reaction[44] and pp — ► p(r]TT + Tr~)p reaction[45] have been compared 
with theoretical works[23, 19, 32, 46]. The former reaction is expected to give data of the KK — > ni] 
process, and the latter to derive information on the elastic nr] — > ttt) process, but the data is not sufficient. 
We do not give the phase shift and cross section here, therefore. (These are shown in Sec. 4.) 

The (KK)i-i scattering amplitude in the single channel formalism has a bound state, but the situation 
is much different from the /o case; the bound state disappears when the irr) coupling is switched off, and 
the irr] elastic channel is repulsive in the single channel formalism. The channel coupling between the irrj 
and KK channels is strong in this channel, however, and then the production rate (1 — J?? 10 ))/4 reaches 
almost the maximum value 0.25. This strong channel coupling generates the ao structure as pointed out 
in the meson-exchange model [31] and unitarizcd chiral theories[32, 33]. Thus, the generating mechanism 
of ao (980) is different from that of / (980); the channel coupling is crucial for both, but it is strong for 
the former but weak for the latter. 




Figure 3: Energy dependence of — Ty. (a) — T\\. (b) — T 12 - (c) — T 22 . The lines are the same as in Fig. 2. 



The amplitudes in the (1,0) channel show the common shape having a sharp cusp-like peak at 
the KK threshold as shown in Fig. 3. This feature will be examined in the -K-q mass spectra observed in 
radiative reactions <fi — > nrj and 77 — ► irrj. 

2.4 Other partial waves 

The (/, J) = (1/2,0) and (3/2) phase shifts are described qualitatively well by our OOP amplitudes, 
though the results are not shown. But better fitting is obtained if we use fx = 1-22/^. The scattering 
lengths arc 

a 1/2) o = 0.216 (0.161) a 3Ao = -0.054 (-0.049) (2.14) 

in units of (m^ 1 ), where the values in the parentheses are for the case that fx — 1.22/^, while ChPT 
gives [47] 

ai/2,0 = 0.17 ±0.02 a 3/2i0 = -0.05 ± 0.02. (2.15) 

The vector resonances, p and if* (890), are also reproduced well by our OOP amplitudes. The mass 
and width of the p meson are 784 MeV and 162 Mev, respectively. Those of the K* meson are 818 MeV 
and 28 MeV, respectively, for fx = fn- In the case that fx = l-22fn the obtained values are 900 MeV 
and 45 MeV, respectively. 

Thus, we say that the validity of the two-channel OOP amplitudes is confirmed. 
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2.5 Pole search for the f and a states 

The OOP amplitudes Tij can be decomposed into a pole term and its background term as 

Ik(*)= T i (s)A(s)T j (s) + U ij (s), 



(2.16) 



where T^s) is the one-particle irreducible form factor of the i-th channel, A(s) is the propagator of the 
relevant particle, and Uij is the background amplitude satisfying unitarity by themselves [48]. The form 
factor Ti satisfies the unitarity relation, 



imr i (s) = -[/*.( s )p j ( s )r j ( s ), 



(2.17) 



and then the phase of Ti is controlled by those of Uij. 

The propagator A(s) continued to the second Riemann sheet has a complex conjugate pair of poles, 
and the pole near the real axis on the lower half-plain contributes to the physical amplitude as 



1 



9je" 



(2.18) 



where gie ai is the form factor T l (s) at s = M R , and Wr = Mr - iMj with Mr (Mj) being the 
real (imaginary) part of the pole. The poles and residues are tabulated in Table 1. 



Jo (980) 


a (980) 


W R = (978.4- 19.7 i) MeV 


W R = (1110.4- 8.9 i) MeV 


.gi = 1.54 (GcV) 


g 2 = 4.57 GcV 


gi = 6.63 GcV 


32 = 8.50 GcV 


ai = -85.6° 


a 2 = -15.8° 


ai = -8.5° 


a 2 = -36.1° 



Table 1: Pole positions, residues and phases of the /o(980) and ao(980) states. 



It should be noticed that the pole of the a state appears above the K K threshold by 120 MeV, while 
the pole of the /o state exists below the threshold and the distance between the two poles is about 130 
MeV. Nevertheless all the amplitudes of the isovector channel show the cusp-like behavior at the KK 
threshold, and then the two peaks appear very close to each other in physical processes. The large phase 
«i for the /o pole indicates a large contribution from the background of the irir — > tttt amplitude, that is 
due to the a enhancement centered at 500 MeV. On the other hand ol\ for the ao pole is small, and then 
n-q — > ttt] scattering may be dominated by the pole term, but the size of a 2 indicate that background 
contribution is not ignored. 

We quote examples of the poles and their residues of previous works: 



fa W R = (0.987 - 0.011 i) GcV (31, g 2 ) = (1.18, 3.83) GeV, 
a W R = (1.030 - 0.086 i) GeV (31, g 2 ) — (4.08, 5.60) GeV 



(2.19) 



in the two-channel calculation [49], and 

fo W R = (0.989 - 0.022 i) GcV (31, g 2 ) = (1.155, 1.227) GeV (2.20) 

in [15], where g 2 is much smaller than that of Eq.(2.19) and our value in Table 2. 
If we define the partial width as an integral[15], 



2Mj_ r Pi (s)gf 



we have 

= 45.3 MeV T kR = 39.4 MeV (2.22) 

for the / (980) state, and the sum of the partial widths, 77.8 MeV, is larger than 2M/ = 35.4 MeV by a 
factor 2. As to the ao(980) state the same definition gives huge partial widths; 

= 567MeV T kR = 554 MeV (2.23) 
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because of the large residues. Such huge partial widths have already been pointed out by Flatte[50], who 
shows that the partial nrj width could be 300 MeV, when the Flatte phase space factor pf{s) is usued, 
where Pf(s) = P2(s) for s > s 2 , and Pf(s) = i\p2(s)\ for s < s 2 . These facts cast doubt not only on the 
above definition of the partial widths but also the conventional resonance interpretation of a in contrast 
to an isolated narrow resonance. 

In order to to extract the pole contribution from the whole amplitudes, we have to express the pole 
term Pjj (s) in a region away from the pole position. After testing the Breit-Wigner form, the Flatte form 
and the loop form used in Ref. [4], we find that the Breit-Wigner form is suitable to the fo state, where 
we compare the Breit-Wigner amplitudes 

Pij{S) = s-M^ + t Sn( P i9 2 1 +P29 2 2 ) (2 - 24) 

with the Tij amplitude; we observe that the two 12 amplitudes resemble with each other near the fo 
resonance peak, except for the phase because of neglecting the phase a\ 2 ~ 90° in P 12 - See Fig. 4. 



al a2 




0.7 0.8 0.9 1 1.1 1.2 0.7 0.8 0.9 1 1.1 1.2 



Figure 4: (a) The fo state. The left side is — P\ 2 and the right side +ReTi 2 , — ImXi 2 . (b) The ao state. 
The left side is —P\ 2 and the right side is — T\ 2 . The solid lines are the absolute values, the dotted lines 
the real parts and the dashed lines the imaginary parts. Note that the real and imaginary parts are 
exchanged between the left and right side of (a) because of a\ + a 2 ~ —90°. 



But it turns out that in order to reproduce the sharp cusp-like behavior in the isovector channel the 
Flatte form is required, where pp is used instead of p 2 , as shown in Fig. 4 (b). The loop form is useful 
to the fo state, but not for ao state, because the depression of the real part is insufficient below the KK 
threshold. The residues strongly depress scattering amplitudes below the KK threshold in the Flatte 
form. 

Another method to extract the pure fo amplitude is to assume the form of Uij as a single channel 
elastic amplitude as 

ti?(*)-ttf(*Wi(*)[tiiW 
and subtract it from Tn. The extracted fo amplitude is given as [51] 



Ms) = (T„(j) - [/ 11 (j))e- M »« (2.26) 
where Sjj is the phase shift of V\\. From this definition the mass and the width of the fo state are 

M fo = 977.7 MeV r /o = 35.6 MeV, (2.27) 
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where Mf is defined as a mass at which the phase shift of fo(s) crosses 90°, and 



T f° ~ 2^/r^= M fo- (2 - 28) 

The values in Eq.(2.27) are very close to the pole position in Table 1. In this case there are no background 
term in tttt — > KK and K K — ► KK amplitudes, and T i2 and T 22 are dominated by the /o state alone. 
Indeed, the clear peak structure appears both in these amplitudes, but it would not be true that there 
are no background term in these amplitudes. 



3 The fo and ao states in production processes 
3.1 Radiative <f> meson decays 

Many theoretical calculations of the radiative 4> meson decays are based on the common production 
mechanism, in which the final S-wave two-meson scattering amplitude is connected to the initial decay 
vertex through the charged kaon loop. The validity of the on-shell factorizability of the loop and the 
scattering amplitudes is discussed in Refs. [32, 6, 7, 52]. 

According to Ref.[7] we write the decay amplitude of the <\> meson to two-meson state jmim 2 as 



F{<f> -> imi m 2 ) = 2eg{G K (s) + f- 



2mj 



-Jk(s)}T k+k - 



(3.1) 



where m\m 2 means 7r°7r°, 7r°?7 or K°K°, and g and / are the parameters defined as 

_ Gym^ j j. F v 



V2.tf 



and / 



2G V 



- 1 



(3.2) 



with Gy and Fy being the constants in the chiral Lagrangian[53, 54, 7] and a factor l/y/2 in g is the 
ratio of the coupling constant pmr to cfiKK. The triangle KK loop integral Gk(s), which connects the 
<f> meson to the S-wave two- meson state after emitting a photon, is given as [3, 55, 5] 



G K {s) = 



mi 



(4tt) 2 



m" A 



s-m" A 



<tk{s) log 



log ; 



ok{s) + 1 
<tk(s) - 1 



1 



log 2 



a K {mMog 



'<rK(mD + 1 



(3.3) 



^K) - 1 



and Jk{s) is the kaon two-point loop function given in Eq.(2.6). The scattering amplitude TV+_ff-_> mim2 
is T K + K -_> nn or T K + K -^ V1} , which is obtained in the previous section. We adopt the values Gy = 55 
MeV and Fy = 165 MeV given in [7], which are suited to the <j> — > K + K~ and (j> — > e + e~ decay widths, 
respectively, and then we have g — 4.69 and / = 0.5. We emphasize that there are left no adjustable 
parameters, and then the mass spectra are a proving ground of the validity of Ti 2 (s) below the KK 
threshold. 

The mass dependence of the two-meson state in the radiative decay is given as 



rfT(s) 
d\fs 



(-) (~ 

\3ttJ V47T 



k f( m l - s ) 



F(<f> — > 7mim 2 ) 



(3.4) 



where kf is the momentum of a final meson in the rest frame of the final {m\m2) state, and we define 
F as F — 2egF. The s-dependence of the phase space factor fc/(m^ — s)/m| has the maximum near 
the middle of the whole y^-range, and vanishes at both ends. Furthermore, the loop integral Gk{s) + 
f( m 4> ~ w2 )/(^ m 1) ' Jk(s) has a strong cusp behavior at the KK threshold, which is very near to m^, 
so that the mass distributions near 1 GeV are affected double by these kinematical factors. 

The calculated mass distributions of the <j> —> 77r°7r° and <f> — > jir i] decays reproduce the experimental 
data fairly well except for the mass region from 400 to 650 MeV of the <j) — > jmr decay, where the mass 
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0.30.40.50.60.70.80.9 1 0.70.750.80.850.90.95 1 



Figure 5: (a) dB{<j> — > 7r°7r°) / 'd\fs x 10 4 /GcV. The solid line is the result of the full amplitude and 
the dotted one that of the pole-terms. Experimental data are taken Experimental data are taken from 
[56, 57, 58]. (b) dB{<j> — > "f^rj)/ 'dyfs x 10 4 /GeV. The lines are the same as in (a). Experimental data 
are taken from [59, 60]. 



distribution shows a shallow dip. See Fig. 5. Our results strongly indicate the validity of both of 
our scattering amplitudes and the production mechanism. The calculated and experimental branching 
fractions are tabulated in Table 2 and 3, respectively. The tttt mass distribution below 700 MeV for 
(p — > 77r°7r° is affected by dominant backgrounds W7r° — > 77r 7r° and pn a — > 77r°7r° as stressed in [56, 57]. 
It is better, therefore, to compare a theoretical calculation with experimental date above 700 MeV. 





B{(j> -> 77T U 7T U ) 


B(4> -> 77r u r?) 


Whole mass range 


1.29 x 10~ 4 


0.875 x 10~ 4 


w > 700 MeV 


0.917 x 10~ 4 




w > 900 MeV 


0.493 x 10~ 4 




Pole contribution 


0.655 x 10~ 4 


0.701 x 10~ 4 



Table 2: Branching fractions integrated over the indicated mass ranges. 



B{4> -» 77r°7r°) • 10 4 


CMD-2 


SND 


KLOE 


whole mass range 


1.08 ±0.17 ±0.09 


1.158 ±0.093 ±0.052 


1.09 ±0.03 ±0.05 


w > 700 MeV 


0.92 ±0.08 ±0.06 


1.034 ± 0.066 ±0.046 


0.96 ±0.02 ±0.04 


w > 900 MeV 


0.57 ± 0.06 ± 0.04 


0.559 ± 0.053 ± 0.025 




5(0 -» 77^77) • 10 4 


CMD-2 


SND 


KLOE 


whole mass range 


0.90 ±0.24 ±0.10 


0.88 ±0.14 ±0.09 


0.851 ±0.051 ± 0.057* 



Table 3: Experimental data on the branching fractions, where CMD-2 is referred to [56], SND to [57, 59] 
and KLOE to [58, 60]. * means the data in which 77 is identified through 77 — > 77, and it becomes 
(0.796 ± 0.060 ± 0.040) when 77 -» 3tt decay is used. 

Our results B{<j) -> 77r°7r°) = 0.917 x 10" 4 above 700 MeV and B(</> -» 7^77) = 0.875 x 10" 4 
are consistent with the data in Table 3. The problem is what parts of the mass distributions can be 
attributed to the contributions from "pure" /o(980) and ao(980) states. If we use the pole terms, the 
usual Breit-Wigner form for /o and the Flatte form for ao given the previous section, we have B{<p — > 
7/0 -> 77T7r) = (0.655 x 3 = 1.965) x 10~ 4 and B(<j) -» 7 a -» 7^77) = 0.701 x 10" 4 , both of which are the 
values integrated over the whole mass range. The /o pole contribution above 900 MeV is 0.413 x 10~ 4 , 
that is 84 % of the 7r°7r° distribution in the same mass range, and B((f> — > 700 — ► 7^77) = 0.701 x 10~ 4 is 
80 % of the whole 7r77 mass distribution. The pole contribution dominates the mass distributions above 
900 MeV for the 77r°7r° decay, while the pole contribution dominates the whole mass distribution for 



9 



the 77r°?7 decay. Here we note that these values can be regarded as the branching fractions of the whole 
B((f> — > 7/0) and B{<j> — > 7a ), because the -RTJ? contributions are negligibly small due to the small phase 
space. 

Finally we obtain 

B{4> -> ~fK a K a ) = 4.26 x 10~ 8 (3.5) 

for the neutral kaon pair decay. In order to get this value the mass difference between the neutral and 
charged kaon plays an important role to suppress the phase space; we use 497.67 MeV for K° here, while 
495 MeV is used for the kaon mass throughout this paper. Our result is very similar to the values of 
previous works; 4.36 x 10~ 8 in the second paper of [4] and 5.0 x 10~ 8 in [6]. The charged kaon pair 
production is dominated by the Born terms including the Brcmsstrahlung terms, and then it is almost 
irrelevant to the /o and do states. 



3.2 Two-photon collision processes 

Two-photon collision processes are not so suitable to study the / states in contrast to the radiative 
(f> meson decays, because the charged kaon loop dominance does not hold. The Born terms with the 
charged pion exchange and u>-meson exchange contribute to 7r + 7r~ and 7r°7r° production, respectively, 
and the charged pion loop and W7r°7r° loop diagrams contribute to these production processes. And there 
are contributions from higher partial waves and the phase space are opened. Indeed, experimental irir 
mass distributions seem not to show manifest signals of the /o state[61, 62, 63]. On the other hand the 
do production seems not to be worried with various production mechanisms besides the charged kaon 
triangle loop diagram as in the radiative <f> meson decay. We discard the W777 and p^ir vertices, because 
their coupling constants are much smaller than that of the u>jir° vertex. Focusing on the S'-wave pion 
pair and ir°r) production processes we study 77 — * m\m2 processes through the charged pion and kaon 
loop diagrams and the cj exchange diagrams, which are developed in Refs.[64, 65, 66]. 

The production amplitude Ff(w) from the initial photon state with the helicity (++) to the final 
S'-wave 77117712 state / can be written 

F f (s) = -2e 2 ^B f (s) + £ T f i{a)Gi{8)j , (3.6) 

where the final state / denotes charged pion pair, neutral pion pair and 7r°?7 states, and the intermediate 
state i = 1 means the charged pair, i = 2 the charged kaon pair state, and i = 3 the 7r°7r° state generated 
in the ujtt^tt triangle loop. Bf(s) is the Born term, but B mi — 0, and T/j is the i —* f scattering 
amplitude. The triangle loop functions G\(w) and G2(s) are the integrals of the charged pion and kaon 
loop diagram, respectively, which is given as 

aw /i + "4vf2j£tm. (3 .-) 



(4tt) 2 [ s \ai(s) - 1 

The loop function G3 with the u> exchange without any form factor is given in Ref . [64] , which is rewritten 
in Appendix because of a lengthy expression. The discontinuity of Gi (s) across the physical cut gives the 
Born term as 

2m 2 {l + <n(sy 



ImGj(s) = -pi(s)Bi(s) = -pi(s) 
for i = 1 and 2, and 



n log . , , , 
sai{s) \l-ai(s) 



(3.8) 



= -*«*<■> - -*<•> ■ ^ {-+ J) '-(Sra)} • <"> 

where B 3 is the u> exchange Born term with = 1.35 GcV -2 [65], v(s) — 2(M 2 — m 2 )/s, and a factor 
l/y/2 is related to the definition of the cross section. The production amplitude Ff satisfies the final 
interaction theorem as 

Imf>( S ) = -Y,Tj i {8)p i {8)F i (8),. (3.10) 
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The S-wave two-meson production cross section of the / state is written as 

2 



with Ff = —2e 2 Ff, and kf is the CM momentum. 



Ff(s) 



(3.11) 




0.40.60.8 1 1.2 




0.40.60.8 1 1.2 




.70.80 



Figure 6: (a) 77 — ► tt + it cross section integrated over | cos9\ < 0.6. The solid line is the total S-wave 
cross section, and the dashed one the Born term. The behavior near the fo region is given in the inner 



box. (b) 77 



cross section integrated over |cos#| < 0.8. The solid line shows the sum of the 



S- wave cross section and the / 2 (1280) resonance. The dotted line shows the 5- wave cross section and the 
dashed one the / 2 (1280) cross section. The experimental data are taken from [61]. (c) 77 — > -K rj cross 
section integrated over | cos#| < 0.9 with the a 2 (1320) resonance. The lines are the same as in (b). The 
experimental data are taken from [68, 69] 



The calculated S'-wave tt + it~ cross section integrated over the region | cos#| < 0.6 shows a large peak 
resulting from the Born term and a small fo signal coming from the rescattering terms as shown in 
Fig. 6 (a). The cross section integrated over the full angle is merely 18.6 nb/GeV at the f peak. This 
is due to destructive interferences among the Born term and various rescattering terms. The small fo 
signal in the charged pion distribution appears to be consistent with the experimental mass distributions, 
which show no clear f signal[62, 63]. The calculated S'-wave 7r°7r° cross section integrated over the region 
I cos 9\ < 0.8 shows a dip and a cusp, but not the fo signal. If we discard both the u> exchange Born term 
and the rescattering term with the ojirir loop, we can see a clear f peak. This implies that the present 
form of the u> exchange amplitudes without any form factor strongly reduces the fo signal. We should 
remember, here, that the Tn amplitude does not show a peak but a sharp dip in the fo region, and that 
the rescattering terms with the charged pion loop and the ivmr loop contain the very Tn amplitude. 
The broad bump around 500 Mev comes mainly from the pion loop diagram with Tn. Our neutral pion 
mass distribution seems not to be inconsistent with the experimental mass distribution, but there may 
be some deficiencies in the region of fo state. We think, however, that these possible deficiencies do not 
imply the invalidity of the OOP amplitudes, but suggest the inadequacy of the treatment of production 
mechanisms. 

As to the 7r°?7 production channel the sharp peak appears in our calculation, but the consistency 
with the experimental data is not clear. It is noted that we take into account the rescattering term with 
the charged kaon loop alone in order to calculated the ir°r) mass distribution. If there contribute other 
production mechanisms, the similar shape of the scattering amplitudes would not destroy the sharp peak 
structure. As stressed in Ref. [66] the axial vector meson exchange would enhance the peak. 

The partial widths of the /o(980) — > 77 and ao(980) — ► 77 have been expected to play a role to 
distinguish the structure of them as in the radiative <j) decays [17]. We have learned in the above that 
different production mechanisms and their interference mask the fo signal in the charged and neutral 
pion pair production. We should be cautious, therefore, of the way how to estimate the partial /077 
width. 

Taking the narrow width limit of the relativistic Breit-Wigner formula given in [61, 62], we can define 
the partial width by integrating the mass spectrum around the peak region as 

T 77 (i? ->/) = _« y dsa^ f (s), (3.12) 
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where ■ s fsl {^Js^) is 0.9 (1.05) GeV for f and 0.8 (1.2) GeV for a . The integration range is chosen so 
as reflect the form of the mass spectra. We obtain 



r 77 (/o -» 7T+TT-) = 0.056 keV, (3.13) 

but we cannot estimate r 77 (/ — » 7r°7r°), because there is no /o signal. If we use the pole terms, Pij(s), 
written by the Breit-Wigner forms instead of , we observe peaks in both spectra deformed by the cusp 
behavior, and we obtain 

r 77 (/ ^7T+7r-) = 0.081 keV, 

r 77 (/o -> ^°7r°) = 0.035 keV. [6 > 

The sum is 0.116 keV, and we have 

T 77 (/o -» all) = 0.20 keV, (3.15) 

if we use B(fo — > 7T7t) = 0.58 calculated from the partial widths Eq.(2.22). 

Though the 7r?7 spectrum cannot be expressed by the Breit-Wigner form, we use the same equation 
with the peak at the KK threshold, 990 MeV, and obtain 

r 77 (a -> TFT?) = 0.17 keV. (3.16) 

If we use (— ImT22 • IG2 (s) | 2 ) instead of (pi(s)|7i2G2(s)| 2 ) in order to estimate the total photo-partial 
width[66, 8], where T22 is T K + K -^ K + K - and T12 is T K + K - ^^a^, we obtain a value, 

r 77 (a -> all) = 0.27 keV. (3.17) 

Due to the similarity of the pole terms to the full amplitudes, the photo-partial widths are almost 
unchanged. 

Preceding theoretical estimates are 



/o(980) \r 77 (/ ^all) = o!28lg;gkeV [67]! (3 ' 18) 

«o(98o) [l^r !!!! : ™™ S (3-19) 



The experimental data using the Breit-Wigner fits are summarized as follows: 

( r 77 (a ^7r?7) = 0.19^'io keV [68], 
a (980) I rTT(a -> tttj) = 0.28 ± 0.04 ± 0.10 keV [69], (3.21) 
[ .r 77 (a -» 7T/7) = 0.24 ± 0.08 keV [70]. 

In the above T 77 (/o — > all) of [62] and [61] are obtained from 7r + 7r~ and 7r°7r° spectrum, respectively, 
using the conventional isospin factor and the old PDG value of B(fo — > 7T7t) = 0.78[71]. It is dubious, 
however, to use a conventional isospin factor to estimate the full nn decay width, because complicated 
interferences should probably mask the /o signal. As to the do signal our value 0.17 keV is not inconsistent 
with the value 0.19 keV of [68], but is a little bit smaller than the average value [70]. 

Finally we briefly comment on J/tp decays. The decays J/tp — > (</> or uj)(-ktt) and J/ip — » p(nri) are 
also expected to specify the quark contents of the /o and ao states [72]. The two-pion mass spectrum of 
the former decay accompanied by <p shows a clear peak at the KK threshold as like as in <p — > 7(7T7t)[73]. 
This strongly suggests that the decay proceeds mainly through the kaon loop so as to make T K ^_, njr 
work effectively. On the other hand the mass spectrum of the two-pion decay with u> shows a large broad 
peak centered at about 450 MeV and a dip-bump structure near the KK threshold [74], that resembles 
to the 77 — ► 7r + 7r~ process. The /o signal seems to be masked. Indeed, the mass spectrum calculated in 
terms of the Born term, that is a direct 5-wave tttt production term, the pion and kaon loop terms and 
a sequential decay mode 6i(1240)7r — ► U7T7T gives the broad peak originated from the Born term and a 
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dip-cusp structure instead of the /o peak^]. 1 Further higher spin meson exchange diagrams would not 
be discarded. 

The J/iji — » p(nri) decay would also include some decay mechanisms similar to the Jji\) — ► w(7T7r) 
decay, but the ao signal might not be masked so much as contrasted with the /o case, since the isovector 
scattering amplitudes have similar shape and do not have large backgrounds. 



4 Large N c behavior of scalar scattering amplitudes 

In the previous two sections we have found that the OOP version of unitarized chiral theories can 
reproduce not only scattering but also production processes fairly well. It is not clear, however, what this 
success tells us about the so-called quark contents of the /o(980) and ao(980) states, because the OOP 
amplitudes contain phenomenological constants, which should be determined so as to be suited to the 
experimental data. 

It is said that ChPT is an effective theory of low energy QCD[76, 25], and qq meson states become 
weakly interacting narrow resonances in the large iV c limit of QCD[77, 78]. It offers a useful means 
to reveal the the quark contents, therefore, to examine whether the scalar mesons survive as narrow 
resonances when N c increases from the physical value N c = 3. (See Ref. [28] for details.) 

1 /2 

Because the pion decay constant f v is of 0(N C ) and the LECs, L\, L2, L3, L5 and L§ are to be of 
0(N C ), but 2Li - L 2 , L 4 , L 6 and L 7 arc of 0(1) [26, 79, 80, 81], we put 

— - N 

L n (N c ) = Ln -^ + AL n , (4.1) 
fl{N c ) = fl~, (4.2) 



where L n satisfy the relations, 2Li — L 2 = L4 = L e = L7 = 0, and AL n are of 0(1). Thus, we have 

L n L, n A_L„ 3 



(4.3) 



ft tt N ° 

with /„■ = 93 MeV. In the case of our LECs AL„ are given as 

AL 2 x 10 3 = -0.10, and AL 7 x 10 3 = -0.25, (4.4) 

and others are zero. The T^{s,t lU ) is of OiN- 1 ) and T^ ly (s,t,u) is also of 0(N~ 1 ), because L n /ff 
scales as 0(N°) as seen in Eq.(4.3). The s-channel loop term given by t ( - 2 \s)J(s)t^(s) is of 0(N~ 2 ). 
This difference of the iV c dependence produces different behavior of the amplitudes when N c becomes 
large. The change of the renormalization scale affects values of AL„ , but we do not consider the scale 
change explicitly because the AL„ terms fade out as N c increases. 

At first, we discuss the behavior of the p meson in the single channel calculation. The mass of the 
p meson is given as a zero of a sum of the term [t^ — t^ ] of OiN^ 1 ) and the real part of the loop 
contribution of 0{N~ 2 ). Since the mass is controlled by a combination of the LECs {2L\ — L2 + L3) 
of the former term [20], the mass stays almost constant for N c > 3. Because the imaginary part of 
the loop term contributes to the width, the width decreases as 0(N C ^ 1 ). Thus, the p meson remains as 
narrower resonance almost at the same point as A^ c increases as shown in Fig. 7. Other vector mesons 
such as if* (890) behave as the same as the p meson. If we extend the calculation to the two-channel 
OOP version, the results do not change. Thus, we can conclude that the vector mesons described by the 
OOP version behave as narrow resonances, and then have the nature consistent with the qq mesons. 



As to the scalar mesons the situation is drastically changed. The N c dependence of the phase shift 
and the cross section of the (I, J) = (0, 0) channel are shown in Fig. 8, where N c increases from 3 to 12. 

1 Recalculation with the present set of the LECS also shows a dip with a rather small cusp. 
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a b 




0.4 0.6 0.8 1 1 .2 0.4 0.6 0.8 1 1 .2 



Figure 7: N c dependence of phase shift (left) and cross section (right) of the p channel. Lines correspond 
to N c — 3, 5 and 15 from the top to the bottom. 

In contrast to the vector channel we observe that the phase shift becomes flat and the cross section fades 
out as N c becomes large; the sharp rise of the phase shift near the K K threshold and the large bump of 
the cross section near 500 MeV seen at iV c — 3 to 5 disappears even at N c = 6, and then the phase shift 
and the cross section become almost fiat and fade out. Similar drastic change in the N c dependence has 
also been observed in Ref. [82], though it is in quite different context. 




0.4 . 6 0.8 1 1.2 0.4 0.6 0.8 1 



Figure 8: N c dependence of the phase shift (left) and the cross section (right) of the (0,0) channel. Solid, 
dotted, and dashed lines correspond to Nc —3, 6 and 12 cases, respectively. 



The (/, J) = (1, 0) channel contains the a (980) state and appears as a cusp-like sharp peak at 
N c — 3 as seen in Fig. 9. The rising phase shift after the cusp bends down and becomes to a flat curve, 
and the cross section having a sharp peak fades out as iV c increases. The (I, J) = (1/2, 0) channel show 
the same behavior as (0, 0) and (1, 0) channels. 



a b 
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Figure 9: N c dependence of the phase shift (left) and the cross section (right) of the (1,0) channel. iV c = 3, 
4 and 12 from the top to bottom. 



We have shown in Sec. 3 that the / (980) pole exists at (978.4 — 19.7 i) MeV and the a pole at 
(1110.4 — 8.9 i) MeV at N c — 3. Where do the poles go away as N c increases ? We make an approximate 
calculation of the pole positions by expanding the amplitudes in powers of ki up to the first order, where 
k 2 is the momentum of the KK channel. We observe that the poles move into the upper half plane of the 
IV sheet from the lower half plane of the II sheet, winding around the branch point at the KK threshold, 
and go away from the real axis as shown in Fig. 10. This approximation gives Wf = (977.9— 18.5 i) MeV 
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and W ao = (1095.7 + 2.7 i) GeV at iV c = 3, which are close to the above values. Pole positions at larger 
N c cannot be reliable owing to the rough approximation, but this behavior would remain unchanged. 



0.2 
0.15 

0.1 
0.05 



-0.05 



Figure 10: N c dependence of the /o(980) pole (left) and ao(980) pole (right). Both of the poles wind 
around the branch point at K K threshold to go upward on the IV sheet. Black point starts at N c = 3 
and increases by 1 to 10 for f and to 6 for a . 



We have calculated the N c dependence of the vector and scalar channels stating from N c — 3 to finite 
values, 30 for the vector channels and 15 for the scalar channels within the two-channel OOP version, 
and we have observed that the vector mesons survive as narrower resonances at almost the same position, 
but the resonant structures of the scalar channels fade out at rather low values of N c near 5 or 6. By 
extending these observations we are led to a conclusion that the light vector meson nonet has the nature 
consistent with the qq mesons in large iV c QCD, but the light scalar meson nonet cannot survive in large 
N c and then cannot have the nature of the qq mesons. This conclusion is the same as obtained by Pelaez, 
excluding an exceptional case of the ao(980) state[27], where the full 0(p 4 ) amplitudes are used. It is also 
consistent with the result that the scalar mesons are describable without pre-existing tree resonance poles 
and then of dynamical origin, possibly except for / (980)[46, 49]. If one treats the /o state within the two- 
channel model, the pre-existing resonance pole is not needed [49]. Our conclusion supports the arguments 
that the light scalar mesons are of the KK molecule[29, 30, 31], and of q 2 q 2 states[83, 84, 85, 3, 4, 5]. 
It is interesting to see that the lattice QCD calculation suggests that the masses of scalar mesons with 
1 = and 1/2 are almost twice of those of the p and if* (890) mesons, respectively, if they are composed 
of connected qq lines [86]. 

5 Concluding remarks 

In this paper we have discussed S'-wave scattering processes and radiative production processes of two- 
pseudoscalar mesons in a unified manner. The scattering amplitudes are constructed within the two- 
channel OOP version and satisfy exact s-channel unitarity. We emphasize that the results of the radiative 
cj) meson decays are the prediction by the two-channel OOP scattering amplitudes, once the decays occur 
dominantly through the charged kaon loop. 

We have found that the scattering data are reproduced rather well, though we do not attempt to get 
the best parameter set. Radiative production processes reveal the structure of the transition amplitudes 
T12 below the KK threshold. Our results of the radiative (f> meson decays show the validity of our 
amplitudes under the dominance of the KK loop mechanism as discussed in Sec. 3.1, where the 7r°7r° 
and tt°t] mass distributions are reproduced fairly well, at least above 700 MeV for the former channel, 
and then the branching fraction B(<fi — > 77r°7r°) above 700 MeV and B((f> — > -yir°ri) are close to the 
experimental data. 

We point out that the branching fraction £?(</> — > 7/0) is strongly model-dependent, though B((j) — > 
77T 7r°) integrated over a mass range given by different experimental groups are mutually consistent 
within errors as shown in Table 3. The SND collaboration assumes that the whole mass distribution is 
dominantly given by the /o(980) resonance with m/ = 969.8 ± 4.5 MeV and r f ^ n7I ~ 200 MeV, and 
gives B{4> -» 7/0) = (3.5 ± 0.3±o j) x 10~ 4 [57]. The CMD-2 collaboration gives B{(j> -» 7/o(980)) = 
(3.05 ±0.25 ±0.72) x 10~ 4 under the single / (980) fit, but (1.5 ±0.5) x 10" 4 under the two-resonance fit 
with / (980) and /o(1200)[56]. On the other hand KLOE collaboration fits the 7r°7r° mass distribution 
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in terms of the j(f + cr(600)) mode with m a = 478 MeV and T a = 324 MeV, which are obtained in [87], 
and gives B(<j> -> 7/0) = 3 x (1.49±0.07) x 1CT 4 = 4.47 x 1CT 4 [58, 60]. In this analysis the a contribution 
gives a bump near 500 MeV, but the experimental mass distribution is very small there, so that the bump 
is erased by the destructive interference with the /o resonance with a large width, where m f — 973 MeV 
and r(/ — ► 7T7r) = 260 MeV. This enlarges the 7/0 branching fraction. It is further pointed out that 
their amplitude violates unitarity[15]. We note that a careful analysis using K-matrix parametrization of 
and D + decays into three pions does not need such a low mass a pole in T-matrix[88]. 

In theoretical estimates performed by adjusting model parameters so as to fit the experimental data, 
the branching fraction over the whole mass range is almost consistent with each other. But theoretical 
estimates of B(<p — > 7/0) seem to be model-dependent similarly to the experimental estimates. Our 
integrated pole contribution B((f> —> 7/0) = 1.97 x 10~ 4 is much larger than the value 0.31 <~ 0.34 x 10~ 4 
estimated through calculating pole residues [15], but smaller than the value 3.11 <~ 3.19 x 10~ 4 obtained 
by integrating ImT 2 2 so as to estimate the /o pole contribution [8]. Thus, the branching fraction of the 
"pure" /o signal has not yet been established in theoretical works originating mainly from uncertainties 
of the parametrization of scattering amplitudes. On the other hand our result B((j) — > ja ) = 0.7 x 10~ 4 
is similar to values of other models. 

Unfortunately, it is difficult to deduce a concrete result on the /o state from two-photon production 
processes, because the processes are not dominated by a single production mechanism, and the large 
charged pion exchange and lo exchange terms mask the /o signal originated from the charged kaon loop. 
Similar complex situation is seen in the J/tp — > ufo process. We have observed that mass spectra 
calculated by the pure pole terms, expressed by appropriate Breit-Wigner forms, arc different from 
those calculated by full scattering amplitudes. This implies that the background coming from the a 
enhancement cannot be discarded in comparing calculated results with experimental data. 

It is impossible to conclude the quark contents of the light scalar mesons even if we succeed in 
describing scattering and production processes simultaneously in terms of the OOP amplitudes. We have 
proposed to study how scattering amplitudes behave as N c increases from 3, because ChPT is an effective 
theory of QCD and has the explicit N c dependence. As a result we have obtained the probable conclusion 
that while the light vector mesons are of typical qq mesons, all of the light scalar mesons below 1 GeV 
cannot survive as narrow resonances when N c increases from 3, and then the light scalar mesons cannot 
be of simple qq states. 

If the mesons in the light scalar nonet are dominantly composed of hadronic or four quark component, 
but include \qq >p with a small fraction as in Ref. [1], we could find out the small \qq >p component by 
increasing iV c in theoretical models, because the large hadronic or four quark component fades out [84] 
and the qq component remains. At least, our calculation within the two-channel OOP approximation 
does not indicate that such an intriguing change will occur in larger N c region. 

Our conclusion strongly indicates that all of the light scalar mesons are dynamical effects originating 
from unitarity, chiral symmetry and channel couplings. 

A Isospin decomposition of scattering amplitudes and a;7r°7r° 
loop integral 

Here we summarize the isospin decomposition of amplitudes written in terms of charged states. 

< 7r + 7T-|T|7r + 7r- > 
< 7r+7r-|T|7r 7r° > 

< 7T 7T |i |7T TT > 

< K + K-\T\ir+ir- > 

< K + K-\T\ir°Tr° > 
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r ii + r ' 

r ll + r ' 
1 

1 

7i 



rriO 
J 1 9 



rpO 
1 127 



<K+K-\T\n°ri> = ±T} 2 , 

<K + K-\T\K°K°> = \{T° 2 -T 2 \). 

These amplitudes satisfy unitarity relations from each other. 

The a>7T 7r° triangle loop function is given in Appendix B of Ref. [64], which is rewritten so as to fit 
our normalization and sign. 

G 3 (a) = (2M^ 2 + S(Tl ( s )log^ l(s) + 1 



\M*-m* \M*-m%, 
L 2 = -Sp{l + iu)-Sp(l-iu)-Sp(-^p-)-Sp{-^—p-) 
~ log(l + ^)log(^)-log(l + ^)log(i-^) 

7T 2 _i 7T , fl+V + ai\ , „ o 

- — +7rtan u + i-log — for s > 2m_, 

3 2 \ 1 + u — <7i/ 

2(M 2 -m 2 ) 
s 



M 2 — m 2 ' 

where Sp(z) is the double polylog function defined as 



Sp(z) = — I dx 
Jo 



log(l - x) 
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